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Abstract — Measurement and modeling of the influence 
function plays a vital role in assessing the performance of the 
continuous MEMS deformable mirror (DM) for adaptive optics 
applications. The influence function is represented in terms 
of Zernike polynomials and shown that the dominant modes 
for representation of central actuators of the DM are different 
from those for the edge actuators. In this paper, a novel and 
effective method of modeling the influence function for all 
the DM actuators is proposed using the 2D sinc-squared 
function. 

Index Terms — deformable mirror, influence function, adaptive 
optics 

I. Introduction 

Deformable mirrors (DMs) are adaptive devices capable 
of correcting the dynamic effects on optical wavefronts due 
to random processes, in the likes of atmospheric turbulence. 
These are often used in adaptive optics systems and are 
commonly controlled by the application of forces on actuators 
that hold the mirror surface. Two kinds of DMs are available 
- segmented and continuous. Unlike the discrete actuators 
holding individual mirror segments in the segmented DM 
case, the later uses a continuous membrane. Minimal edge 
diffraction in the case of a continuous DM makes it more 
desirable of the two. Since in most applications, the 
wavefronts to be corrected by the DM are often much 
smoother than the resolution limit of the device, it is 
advantageous to use the continuous DM with precise 
calibration [1-3]. The effect of individual actuators is desired 
to be localized, but, in the case of a continuous DM, the 
application of a certain force on one of the actuators results 
in a finite displacement (that cannot be neglected) at the 
neighboring actuator location. This error due to non-localized 
actuation can be largely compensated by determining the 
influence function, which is a measure of the response of the 
DM surface when a known voltage (called poke voltage) is 
applied to a single DM actuator (a constant bias voltage is 
applied to all the other actuators). The motive of the present 
work is to construct a mathematical model for the measured 
influence function. 

Two other metrics used for the assessment of the DM 
performance include, the inter-actuator coupling coefficient 
and the interaction matrix. The inter-actuator coupling 
coefficient is defined as the ratio of the displacement of the 
DM surface at the poked actuator to that obtained at its 
neighboring actuator. It is often practically convenient to 
employ the interaction matrix, that relates the wavefront sensor 
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measurements and the weighting associated with the basis 
functions used to fit the wavefront shape on the DM, to 
account for the actuator coupling [4]. 

Obtaining a fit function for the measured influence 
function for all the actuators individually can be useful in 
accurately modeling the DM. Most researchers use a Gaussian 
model, 

IF(x,y) = a + b exp{-(w x x 2 +w y y 2 )} (1) 

which is not a true representation of the influence function. 
The outer portion of the influence function is different from a 
Gaussian as was pointed out by Wenhan Jiang, et. al. [5]. 
With azimuthal and radial modification to the Gaussian 
function, the influence function can be modeled more 
accurately [6]. But, it should be noted that the edge actuators 
are poorly modeled by the earlier suggested models. J. H. 
Lee, et. al [7] used the ID sine function to model the ID 
influence function. In this paper, it is shown that the 2D sinc- 
squared function and Zernike moments can be used to 
effectively model the influence function. The measurements 
were made using the Shack Hartmann Wavefront Sensor 
(SHWS) for all the actuators. The SHWS measures the local 
slopes of the incident wavefront (in terms of the shift in the 
focal spots of the microlenses from ideal positions) and the 
shape of the wavefront is calculated using the least square 
techniques. 

In our study, the Multi-DM, manufactured by Boston 
Micromachine Corporation (BMC) was used. The following 
section discusses the method used for the measurement of 
the influence function. Evaluation of the inter-actuator 
coupling is presented in section III. Section IV details the 
modeling of the influence function. 

II. Influence Function Measurement 

A. Experimental details 

The optical arrangement for the measurement of the DM 
influence function is shown in Fig. 1 . A 15mW Melles Griot 
He-Ne laser was used as the source of light. The laser beam 
was filtered using the spatial filter (S.F) and then collimated 
with a doublet lens (CL). The beam was allowed to reflect 
from the DM surface. The beam was then filtered using the 4f 
lens geometry (L , L 2 are equi-focal lenses, f= 12.5cm). The 
assembled SHWS was made-up of a contiguous microlens 
array (from Flexible Optical BV, APO-Q-P200-F40 at 633nm, 
with a focal length of f=40mm and a pitch of p=200um) and a 
compact, high-resolution progressive monochrome CCD 
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camera (Pulnix TM-1325CL). The camera was placed at the 
focal plane of the microlens array (a distance f=4cm away). It 
is important to note here that the performance of the SHWS 
is highly sensitive to placing the CCD at the focal plane of 
the microlenses [8]. 




CCD 



Figure 1. Optical arrangement for the measurement of the DM 
influence function 

In view of the fact that the clear aperture on the DM is 
4.9x4.9 mm 2 and the pitch of the microlens array is 200 um, 
25x25 microlenses were used to incorporate the entire active 
region of the DM, for sensing. Also, the pitch of the CCD 
pixels is 6.45um, which implies that a single mircolens 
corresponds to 3 1 pixels on the CCD, given the pitch of the 
microlens array to be 200um. For the reason that, the light 
levels and frame rates are not the points of concern, the high 
readout noise (50 e per frame per pixel) is not an issue in our 
case [9]. In contrary, the high resolution of the camera 
provides high-quality, accurate wavefront sensing 
capabilities. The maximum voltage that was applied to the 
DM actuators was limited to 40V to avoid SHWS cross-talk. 
The 140 actuators of the DM are numbered as shown in Fig. 
2 (when quoted, they are written in bold throughout the 
paper). In order to differentiate central actuators from the 
ones at the edges, they are shaded in blue. The corners (1, 
12, 133, 144), although numbered, are positions of zero 
activity and no actuator is present here (shaded black). The 
actuator positions shaded in green are called the penultimate 
actuators here. The behavior of these actuators will be soon 
shown to be somewhat different from center and edge 
actuators. 
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Figure 2. Actuator Map: Blue region is the central portion and the 
blackened actuators are the corner actuators, the uncolored region 
shows the edge actuators and green shaded region depicts 
penultimate actuators 



B. Influence function calculation from SHWS slopes 

The influence function was calculated from the SFfWS 
measurements, assuming the slope sensing geometry of 
Southwell [10] and applying the least square technique. The 
measured influence functions for the edge, penultimate and 
central actuators are shown in Fig. 3. It can be observed that 
the influence functions for each type of actuator (central, 
edge and penultimate) are different from one other. This is 
attributed to the presence of the edge and penultimate 
actuators near the fixed periphery. The systematic error that 
can be observed on the background in Figs. 3 a, 3b and 3c 
can be remarkably reduced by taking the slope values, 
measured when OV is applied to all the actuators, as reference. 
The background-error-free influence function measurements 
are shown in Figs. 3d, 3e and 3f. 




(.1) (b) (c) 
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Figure 3. Measured influence function: A poke voltage of 30V is 
applied to actuators: (a, d) 7, edge actuator (b, e) 19, penultimate 

actuator (c, f) 67, center actuator, with zeros bias voltage. The 
surfaces, complete wavefront reconstructions of the DM surface, 
shown above are each normalized from zero to unity 

It is to be remembered that there can be multiple errors 
involved while sensing the displacement of the DM surface. 
Apart from the experimental errors due to beam alignment 
and control, the only other error arises from the inaccuracies 
in the wavefront sensing algorithm. The experimental errors 
can be largely minimized by adopting cautious beam control 
steps. The wavefront sensing algorithm is limited by the 
wavefront sampling geometry, centroiding algorithm used to 
detect the position of the SFfWS spots and the method used 
for reconstructing the DM surface from slope measurements. 

III. Evaluation of the Inter- Actuator Coupling 

The inter-actuator coupling for all the central actuators 
was evaluated. A histogram of the inter-actuator coupling 
coefficient (computed for the nearest neighbor) at different 
voltages for central actuators is shown in Fig. 4. It can be 
noted here that with increasing poke voltage, the inter- 
actuator coupling coefficient increases. The inter-actuator 
coupling coefficient for the next nearest neighbors is very 
low in the case of central actuators (see Fig. 5). The measured 
influence function shown in Fig. 3 hints a possible asymmetry 
in the inter-actuator coupling which is pointed out more 
convincingly in Fig. 6 for the case of central actuators. Clearly, 
the actuator coupling is stronger in 'x' direction. 
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Figure 4. Inter-actuator coupling for the nearest neighbor at 
different poke voltage. X-axis shows the inter-actuator coupling 
and Y-axis shows the number of occurrences, analysis on central 
actuators alone 




0,2 0.4 0,2 0.4 

Figure 5. Inter-actuator coupling coefficient for the next nearest 
neighbor at different poke voltage. X-axis shows the inter-actuator 
coupling and Y-axis shows the number of occurrences, analysis on 



central actuators alone 
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Figure 6. Histogram showing the asymmetry in the inter-actuator 
coupling along x (marked left, right) and y (marked top, bottom) 



directions 
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IV. Modeling the measured Influence Function 

A. Fitting using Zernike polynomials 

The measured influence function can also be exactly fitted 
using the Zernike polynomials, Z(x,y) via the calculation of 
Zernike moments, a''s: 

IF{x,y) = Y j a i Z i {x,y) ( 4 ) 

Since it is not practical to calculate infinite moments, the 
influence function is approximated using 'AT number of finite 
Zernike moments. The error involved in representing the 
influence function using a finite number of Zernike moments 
can be made sufficiently low by choosing a large enough 
value for 'M' . Influence function for actuator 66 represented 
using M=15 Zernike moments is shown in Fig. 7. The residual 
error (see Fig. 7c) is evaluated by subtracting the measured 
influence function (Fig. 7a) from the fitted model (Fig. 7b) for 
the influence function. 




<■) (b) (C) 

Figure 7. Fit result of actuator 66 with 75 Zernike polynomials 
with application of a poke voltage, V p =20V, zero bias (a) Measured 
influence function (b) Fit function (c) Residual error 

Zernike polynomials are defined on the unit disk and the 
measured influence function is spread over a square grid. 
This shape incompatibility compelled us to calculate the root 
mean squared error (RMSE), that determines the fitting 
accuracy, over one half of the square grid, close to the center 
alone, excluding the circular boundary where the error is 
exceedingly large. The value, log 10 (RMSE) reduces almost 
linearly with the number of computed Zernike moments up to 
j ~ 40 as can be seen in Fig. 8. Beyond j > 40, the reduction 
continues, but at a smaller rate. This is because the dominant 
aberrations in the influence function are the lower order 
Zernikes. 

A plot of the Zernike moments for actuators 55, 67 and 
77 is shown in Fig. 9. Here, it must be observed that Zernike 
moments with index, j = 4, 12 and 24 dominate over other 
Zernike moments for the central actuators. This is because 
the measured influence functions for central actuators are 
circularly symmetric, and so are the Zernike polynomials, j = 
4, 12 and 24 (see Fig. 12). In the higher spatial frequency 
range, j = 40, 60 appear to be more different from 0, than the 
remaining higher order moments. The major drawback in 
representing the influence function with Zernike polynomials 
is that the wavefront sensing is performed on a square grid 
and the representation is done on a unit disk, as is the 
property of Zernike polynomials. This becomes even more 
significant while trying to represent an edge actuator. 
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Figure 8. RMSE reduces with increasing number of Zernike 
moments used for fitting: case of 20V applied to actuator 66 
(central actuator), zero bias voltage 
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Figure 9. Zernike coefficients on fitting IF for actuators 55, 67, 79 
with a poke voltage of 25V, at zero bias 

The edge actuator can be represented with very low residual 
error using the Zernike moments as shown in Fig. 10. The 
dominating moments in the case of edge actuators are j =1,5 
and 6 (see Fig. 11, 12). 

a on 



Figure 10. Fit result of edge actuator, 7 with Zernike moments: case 
of poke voltage of 25V, at zero bias. The residual error is nearly 
zero in the region where the fitting is done 
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Figure 11. Zernike coefficients on fitting IF for actuators 6, 7, 
with a poke voltage of 25V, at zero bias 
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Figure. 12. Zernike polynomials that play vital role in IF 
representation 

B. Sinc-squared fit 

From the observation of the line profile taken along the 
X-axis of the measured influence function (IF) for actuator 
79 as shown in Fig. 13, it can be concluded that the influence 
function can be modeled as a sinc-squared function. This 
novel influence function model can be very effectively used 
to precisely represent the response of individual actuators. 
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Figure 13. Line profile (ID view) of the measured influence 
function 

In the 2D case, a most general sinc-squared model of the 
influence function model can be expressed, 



i T7 model / \ 

It (x, y) = a + a. 



sine 
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or by neglecting the cross term, it can become, 



IF" 



1 (x, y) = a + a t 



sine 
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(3) 



where, a , a ; , s_, c , s 
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Cy s , and c/ are the fitting parameters. 
Looking at the Eqs. (2) and (3), the assumed model for the 
influence function certainly looks complicated with the need 
to estimate multiple parameters, nine and seven in number 
respectively. The expression in Eq. (3) can be further simplified 
by considering c and c to be identically equal to zero, in 
those cases where this approximation to the influence 
function remains valid. The trust-region reflective Newton 
method (large-scale) was used for the optimization of the 
fitting parameters. A central actuator, 78 was fitted using this 
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method and it was found that the residual error is very low 
(see Fig. 14). R-square, that measures the success of the fit 
via the calculation of the correlation between the measured 
influence function and the fitted model, is used as a metric 
that determines the fit accuracy and in this case, it is 0.9901 . 
The motivation behind favoring sinc-squared instead of sine 
function is that with the former, the fitting root mean square 
error (RMSE) is lesser by a factor of 5, than that of the later. It 
was also noted that, in comparison with the Gaussian fitting, 
[see Eq. (1)] the sinc-squared fitting [see Eq. (2)] improves 
the R-square value by a factor of 1.8. In the case of central 
actuators, the cross term in Eq. (2) can further improve the 
fitting accuracy, but the increase in fit accuracy is too small 
(R-square is greater by 0.02) and does not significantly affect 
our modeling and hence may be ignored and Eq. (3) can be 
used for simplicity. 

The values of 'x' and 'y' were chosen to be positive 
integers alone, while using the sinc-squared model of 
influence function, in the case of edge actuators (see Fig. 
15). It was observed that retaining the cross term in the case 
of edge actuators and penultimate actuators leads to 
erroneous fitting and hence Eq. (3) was used for edge actuator 
modeling. 
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Figure 14. Fit result of central actuator 78 with sine functions with 
application of 20V, at zero bias 




Figure 15. Fit result of edge actuator 7 with Eq. (3): case of poke 
voltage of 25V, at zero bias. Here, c , c take non-zero values 

It was found that for the central actuators, the median of the 
R-square value for fitting the sin-squared function in Eq. (2) 
is just above 0.98. Also, for the central actuators, the medians 
of a , c x , c, c m , are 0.0388, 0.0145, 0.0276 and -0.0049 and 
hence may be approximated to zero. This observation allows 
us to reduce Eq. (2) into a 5-parameter optimization problem. 
From our observations, the most probable model of the DM 
using the sinc-squared function as defined in Eq. (2) is: 



IF m " ic '(x, y) = 0.63 



sincl [ x sine - 

0.74 J { 0.92 



x sine] — ^— | + 0.27 
0.65, 



Conclusions 

In this paper, a novel model for the measured influence 
function of a continuous DM is presented. The sinc-squared 
function can be used to closely model the influence function. 
It is indicated that the behavior of the edge actuators and 
penultimate actuators is not the same as that of the central 
actuator which is made clear with the domination of different 
Zernike modes in each of the cases. The inter-actuator 
coupling coefficient increases as the poke voltage increases. 
The most affected actuator is the nearest neighbor and the 
next nearest neighbor is affected a little too. There is 
asymmetry in the inter-actuator coupling coefficient in 
different directions. Actuator modeling is better with sinc- 
squared function due to the requirement of circular aperture 
for representation in terms of Zernike polynomials. The 
current analysis will gain more significance when the DM is 
used in an optical layout for wavefront correction 
accommodating the modeled influence function. 
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